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Heavy Nuclei, from RHIC to the Cosmos 

Spencer R. Klein®' 

^Lawrence Berkeley National Laboratory, Berkeley, CA, 94720, USA 

Ultra-relativistic heavy ion collisions produce a high-temperature, thermalized system that may mimic the 
conditions present shortly after the big bang. This writeup will given an overview of early results from the 
Relativistic Heavy Ion Collider (RHIC), and discuss what we have learned about hot, strongly interacting nuclear 
systems. The thermal and chemical composition of the system will be discussed, along with observables that are 
sensitive to the early evolution of the system. I will also discuss the implications of the RHIC results for cosmic 
ray air showers. 


1. Introduction 

The Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven National Laboratory (BNL) collides 
ultra-relativistic ions at energies up to 200 GeV 
per nucleon. The nucleon-nucleon reactions are 
energetic enough that perturbative QCD is ex¬ 
pected to be able to describe much of the collision 
dynamics. 

The goal of ultra-relativistic heavy ion colli¬ 
sions is to study the properties of matter at ex¬ 
tremely high temperatures and/or densities, with 
an eye to mimicing the conditions present in the 
very early universe, lO/iS after the big bang. 
A specific goal is to look for the Quark-Gluon 
Plasma, a state of matter whereby the protons 
and neutrons in a nucleus ’dissolve’, producing a 
gas of free quarks and gluons. 

These collisions may also be similar to those 
produced when heavy-ion cosmic rays hit the at¬ 
mosphere. In the target frame, RHIC projectile 
gold nuclei have a total energy of 4.3 PeV (20 TeV 
per nucleon), energetically reaching the knee of 
the cosmic ray spectrum. 

Relativistic heavy ion collisions were initially 
studied at the Berkeley Bevatron, SIS and the 
Dubna Nuclotron|Q. More recently, there have 
been higher energy studies at the BNL AGS and 
the CERN SPS. The SPS data is often used as a 
lower-energy comparison point; the SPS collided 
lead on lead, at a center of mass energy of 17 GeV 
per nucleon. The earlier studies found several in¬ 


teresting phenomena. These include: 

Anisotropic flow: Heavy ion collisions may be 
described at least partly in terms of fluid dynam¬ 
ics; the system shows fluidlike behavior [^. 

Strangeness enhancement: Production of 

strange particles is several times larger than 
would be expected from superimposed pp 
collisions!^ [Q . 

J/^/ suppression: Production of J /ip particles is 
suppressed compared to the production of Drell- 
Yan dileptons . 

The latter two observations have been pro¬ 
posed as signatures of the Quark Gluon Plasma. 
However, both phenomena might be due to 
normal hadronic interactions, with the addi¬ 
tional strangeness produced in secondary reac¬ 
tions among the produced hadrons, and the J /p 
suppression due to interactions with the initial 
state nucleons and the other hadrons produced 
in the collision. 

In low energy heavy ion collisions, the inter¬ 
acting baryons stop when the nuclei collide. As 
the collision energy increases, the nuclei gradually 
become transparent, the baryons retain some of 
their initial momentum, and the net baryon den¬ 
sity of the produced system drops. At RHIC, the 
net baryon density at mid-rapidity (near the sys¬ 
tem center of mass) should be near zero. 

This writeup will discuss heavy ion collisions at 
RHIC, starting with observables that probe ther¬ 
mal freezeout, such as the global event charac¬ 
teristics, system size, particle spectra, and non- 
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isotropic flow. Next, the composition at chemical 
freezeout will be discussed, followed by signatures 
of the early evolution, focusing on high pT parti¬ 
cles and charm production. 

2. RHIC 

RHIC is a 3.8 km circumference double-ring ac¬ 
celerator which can collide gold ions at center of 
mass energies of up to 200 GeV/nucleon at a lumi¬ 
nosity of 2 X 10^®/cm^/s, corresponding to about 
1,500 hadronic collisions/sec. RHIC can also ac¬ 
celerate lighter ions. The maximum energy per 
nucleon depends on the charge to mass ratio; for 
protons, the maximum center of mass energy is 
500 GeV. It also collides polarized protons, to 
study the spin structure of the nucleon. The lu¬ 
minosity depends on the species; for protons the 
luminosity can reach 1.4 x 10^^/cm^/s, or about 
700,000 hadronic interactions/sec. 

In the year 2000, RHIC collided gold nuclei at 
an energy of 130 GeV/nucleon. Most of the re¬ 
sults presented here are from this run. In 2001/2, 
RHIC collided gold nuclei at 200 GeV/nucleon, 
briefly reaching the design luminosity, and col¬ 
lided polarized protons, with up to 25% polariza¬ 
tion. The long term program will include stud¬ 
ies with lighter ions, gold-gold collisions at lower 
energies, and deuterium-gold and/or proton-gold 
collisions. 

RHIC is instrumented with two large detec¬ 
tors, STAR and PHENIX, and two smaller exper¬ 
iments, BRAHMS and PHOBOS. A third small 
experiment, pp2pp, studies proton-proton elastic 
scattering!^. The collaborations have very differ¬ 
ent strategies for studying ion collisions. 

PHENIX is designed to look for relatively rare 
observables that are sensitive to the early phases 
of the collision, such as charmed hadrons, J/'0 
and direct photons ||]. The detectors are opti¬ 
mized for particle identification, especially lep¬ 
tons and photons. PHENIX has a two-armed cen¬ 
tral spectrometer; each arm is instrumented with 
charged particle tracking, time-of-flight (TOE), 
a ring imaging Cherenkov counter (RICH), and 
electromagnetic calorimetry. Each arm covers 
a solid angle of 135 degrees in azimuth by 0.3 
in pseudorapidity, where the pseudorapity ij = 


— ln[tan(0/2)], with 9 the particle angle with re¬ 
spect to the beam axis. The center of mass is at 
r] = 0. Forward and backward muon detectors 
cover 1.2 < |? 7 | < 2.2 for muons with momentum 
p > 2 GeV/c. Specialized triggers and a high 
rate DAQ system will collect large samples of the 
selected rare probes. 

The Solenoidal Tracker at RHIC (STAR) is op¬ 
timized to study hadrons over a very large solid 
angle, including multi-particle correlations, and 
measure global event characteristics!^. 

STAR tracks charged particles with r]\ < 1.5 
in a large time projection chamber (TPC) in a 
5 kG solenoidal magnetic field. A silicon vertex 
detector covering !?;! < 1 and two forward TPCs 
covering 2.5 < \r]\ < 4.0 complete the tracking 
system. Strange particles like Ks, A, 5 and ft 
are detected by reconstructing secondary vertices. 
Energy loss in the TPCs and SVT and small TOE 
and RICH systems provide particle identification, 
along with an electromagnetic calorimeter. STAR 
records a great deal of information on each event, 
but can only record data from selected events at 
rate slower than PHENIX. 

PHOBOS records charged and neutral particle 
production over most of phase space, up toby! < 
5.4, to search for anomalous event shapesQ. It 
has two small charged particle spectrometers with 
TOE systems for particle identification. 

BRAHMS is composed of precision central and 
forward spectrometers with tracking and parti¬ 
cle identification, along with counters to measure 
charged multiplicity!^. 

The 4 experiments include identical zero degree 
calorimeters (ZDCs) to measure forward neutrons 
from nuclear fragmentation. The ZDCs are 
intended to provide a common method for lu¬ 
minosity and centrality (impact parameter) mea¬ 
surements, in order to facilitate comparisons be¬ 
tween the four experiments. 

3. Ultra-Peripheral Collisions 

Before discussing central collisions, it is in¬ 
teresting to consider ultra-peripheral collisions 
(UPCs), interactions at large impact parameters 
b (minimum ion-ion separation) where only pho- 
tonuclear and two-photon interactions are possi- 
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Figure 1. pr spectrum of 2 track events observed 
at 130 GeV in STAR. The peak at low pt is char¬ 
acteristic of coherent coupling to both nuclei, as 
expected for coherent photoproduction [^ . 


ble. UPCs can probe a wide variety of physics |p^, 
ranging from electrodynamics in very strong elec¬ 
tromagnetic fields to meson spectroscopy to mea¬ 
surements of gluon shadowing in heavy nuclei to 
tests of quantum mechanics [p^. Photoproduc¬ 
tion of heavy hadrons and quarkonium is 
sensitive to the gluon density in the nucleus, and 
hence to gluon shadowing. 

At RHIC, mutual nuclear excitation (includ¬ 
ing both Coulombic and hadronic interactions) is 
used as a luminosity monitor]!^. The process has 
a large cross section, about 11 barns, and small 
backgrounds. 

One easily observable UPC process is exclu¬ 
sive coherent photoproduction of vector mesons, 
Au + Au Au -I- Au + . These events are char¬ 
acterized by an almost empty detector, contain¬ 
ing only two tracks, with a total event transverse 
momentum px < 2fi/ « 100 MeV/c. The low 
Pt, characteristic of the coherent photon emis¬ 
sion and scattering, is a distinctive signature, as 
data from STAR shows in Fig. |^. At 130 GeV, 
STAR measures a(Au + Au Au + Au + p'^) = 
460 ± 220 ±110 m b[]^ , in agreement with theo¬ 
retical predictions |19||. The ratio of to direct 
7r'''7r“ production is consistent with that mea¬ 
sured in 7 p interactions. 



Figure 2. Schematic view of a heavy ion interac¬ 
tion, showing the different stages of the reaction. 


4. Hadronic Collisions 

Hadronic collisions occur in several stages, as 
is shown in Fig. The nucleons collide, and 
their partons interact. The produced particles 
interact and form hadrons (hadronize). As the 
interactions continue and the number of particles 
grows, the system expands and cools. When the 
average particle energy is low enough, inelastic 
hadron production stops, a transition known as 
chemical freezeout. Slightly later, the interpar¬ 
ticle separation is large enough that even elastic 
interactions cease; this is thermal freezeout. 

The key question in this picture is whether 
the produced particles interact as hadrons (i.e. 
a hadron gas) or as partons (i.e., a quark-gluon 
plasma). Partons produced in the initial interac¬ 
tions may remain free for long enough to interact 
with each other and equilibrate, forming a quark- 
gluon plasma. Or, they might immediately form 
hadrons (hadronize), and the interacting system 
will be a hadron gas. Or, they could initial in¬ 
teract as a quark-gluon plasma, and, then, as the 
system cools, hadronize to form a hadron gas. 

Most of the theoretical guidance regarding the 
quark-gluon plasma comes from lattice gauge the¬ 
ory (LGT). Recent LGT calculations indicate 
that the phase transition between hadron gas 
and quark-gluon plasma, if it occurs, is weak (at 
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Figure 3. (a)Relationship between forward neu¬ 
trons (ZDC energy) and charged multiplicity (the 
charge in the beam-beam counters, Qbbc), as 
measured by the PHENIX collaboration at 130 
GeV for 4 impact parameter bins, (b) The over¬ 
all charged particle multiplicity, da/dNch (solid 
dots), with calculations of the multiplicity dis¬ 
tribution for the same 4 impact parameter bins. 
From Ref. pl|] . 


least second order), and occurs at a temperature 
of 150-200 MeV and an energy density Cc « 1 
GeV/fm^ [p0[ . This calculation is for an infinite 
medium with an infinite lifetime; edge effects and 
formation time are not considered. Although the 
expected system lifetime is only w 10“^^s, calcu¬ 
lations indicate that equilibration occurs quickly, 
so a clear phase change is possible. 

Although it is a key parameter in heavy ion 
collisions, the impact parameter b is not directly 
observable. We use two classes of observables to 
infer the impact parameter. The first is the num¬ 
ber of forward (zero-degree) neutrons. These neu¬ 
trons come from the non-interacting part of the 
nucleus. Enough energy propagates from the ini¬ 
tial collision to dissociate the non-interacting part 
of the nuclei into neutrons, protons and small nu¬ 


clear fragments. The other observables are sensi¬ 
tive to the number of interacting nucleons. Exam¬ 
ples are the charged particle multiplicity or trans¬ 
verse energy. A model is necessary to relate these 
observables to the impact parameter. To avoid 
systematic uncertainties, events are often sorted 
by centrality (i.e. by charged multiplicity), and 
divided into classes, such as the 10% most central 
(those with the smallest impact parameter). 

Figure || shows the relationship between the 
number of forward neutrons (measured in the 
ZDGs) and the charged multiplicity!^. The 
charged multiplicity rises continually as the im¬ 
pact parameter decreases. However, the number 
of forward neutrons is largest at moderate impact 
parameters. In very central collisions, most of the 
nucleus interacts, leaving few remnant neutrons, 
while in very peripheral collisions, some of the 
nucleus remains intact, reducing the number of 
forward neutrons, producing the curve in Fig. |^. 

5. Thermal Fheezeout 

The particles present at thermal freezeout are 
those observed in the RHIC detectors, and are 
relevant for comparison with models of heavy ion 
collisions. The charged particle multiplicity is 
shown as a function of pseudorapidity rj in Fig. 
The multiplicity dN/dr] is roughly flat for \rj\ < 2. 
This central plateau shows that there is boost in¬ 
variance. Within this region, the system appears 
invariant with respect to the longitudinal boost 
(velocity); the expansion may be treated in 2 di¬ 
mensions. 

At 130 GeV, the maximum dN/di] is about 
570, rising to 650 at 200 GeV. This corresponds 
to total multiplicities of about 4100 ± 210 and 
4960 ± 250 respectively!^. These multiplici¬ 
ties are considerably lower than most pre-RHIC 
predictions!^, and seem to be best fit by mod¬ 
els based on a combination of hard interactions 
(calculated by perturbative QCD) and soft inter¬ 
actions (extrapolated from lower energies). Most 
popular cosmic ray air shower codes predict con¬ 
siderably larger multiplicities . 

The dN/dy per participant (nucleon involved in 
the collision) are about 40% higher than at lower 
energies, and also 50% higher than in pp collisions 







5 



Figure 4. Charged multiplicity dN/drj for differ¬ 
ent centrality bins at 130 GeV. A flat plateau 
is visible at mid-rapidity. From the PHOBOS 
collaborationH^ . 


at comparable energiespSf. The multiplicity per 
participant rises smoothly as the number of par¬ 
ticipants increases. 

The PHENIX collaboration measured a trans¬ 
verse energy dEx/dij « 578 GeV per unit rj for 
the 2% most central collisions. This is the en¬ 
ergy released in the collision. In a boost-invariant 
picture (supported by the existence of a central 
plateau), particles emitted into a pseudorapidity 
region St] cine from a region of longitudinal size 
SrjcT, where r Ri h/pQ is the time required for the 
initial interactions to occur. Here, po is the en¬ 
ergy scale for the initial particle production. The 
initial energy density depends on this scale. We 
will use a conservative po = 200 MeV/c, so r 1 
fm/c. The initial volume is irR^SrjT and the en¬ 
ergy density e is 


1 dEj' 
ttR^t dr] 


4.5 GeV/fm^. 


( 1 ) 


This is much larger than the 1 GeV/cm^ that lat¬ 
tice gauge calculations predict is required to form 
a QGP. Although the time scale is not completely 
fixed, it seems hard to stretch t enough to reduce 
e below 1 GeV/cm^. 


The baryon:antibaryon ratio at freezeout is also 
of interest. Some baryons are initially present in 
the gold nuclei, while the rest are produced via 
baryon-antibaryon pair production. Antibaryons 
come only from the latter source. At 130 GeV, 
the p : p ratio is 0.6 ± 0.02 ± 0.06, rising to 0.73 ± 
0.03 for A : A and 0.82 ± 0.08 for S : 5 where the 
first (usually only) error is statistical]^. 

These ratios are quite close to 1. The cen¬ 
tral region is nearly baryon free. Pair produced 
baryons outnumber initial state baryons by more 
than 2:1. However, the net baryon number is 
not zero, showing that there is indeed substan¬ 
tial baryon stopping; many initial state baryons 
are transported over 6 units of rapidity. 

The size of the system at thermal freezeout 
has been measured with 2-particle interferometry 
(Hanbury-Brown Twiss interferometry), taking 
advantage of the Bose statistics that increase the 
abundance of particle pairs with momentum dif¬ 
ference Ap = Pi —p 2 < h/R. Here R is the radius 
of the last elastic interaction. The momentum 
difference vector is decomposed into longitudinal 
(along the beam direction), side (transverse to 
the observer, and out (toward the observer) com¬ 
ponents (the Bertsch-Pratt decomposition). The 
source radius for a source with an assumed Gaus¬ 
sian is about 6 fermi in all 3 dimensions]^. This 
is about twice the initial nuclear radius (about 
6.5 fermi, in a Woods-Saxon [almost hard sphere] 
density distribution). This source size is similar 
to that observed in much lower energy collisions 
at the SPS; the lack of growth is a surprise. It is 
also surprising that the source radii in all 3 dimen¬ 
sions are similar; this indicates that the particles 
are emitted in a very short time scale, i.e. that 
thermal freezeout occurs quite suddenly over the 
entire nucleus. 

Finally, we can compare the charged pion, pro¬ 
ton and kaon spectra. For pr < 2 GeV/c (the 
non-perturbative region), all three spectra are 
consistent with thermal emission. However, the 
three species have rather different apparent tem¬ 
peratures. The different temperatures may be 
due to a collective outward motion known as 
radial flow. If the expanding particles interact 
strongly, they tend to move outward at the same 
velocities. Then, the thermal fit temperature Tapp 
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Figure 5. Temperatures T and expansion velocity 
fi observed at different collision energies. From 
Ref. |§. 


for a particle with mass m is 

Tapp = T + m/3^ (2) 

where T is the actual temperature and /3 • c is the 
collective expansion velocity. The 3 species sat¬ 
isfy Eq. I for T = 120 MeV and (3 = 0.52 [||. 
The system expands outward at more than half 
the speed of light! As Fig. || shows, the temper¬ 
ature is slightly lower than that measured at the 
SPS, but (3 is significantly higher. The tempera¬ 
ture is comparable to the transition temperature 
predicted by LGT calculations. The very large [3 
is characteristic of explosive expansion, with very 
high pressures and strong rescattering. 

Anisotropic flow is another observable. In a 
non-central collision, the reaction zone is ellip¬ 
tical. Pressure converts this spatial asymme¬ 
try into a particle density/momentum anisotropy 
which is usually parameterized as 

dN 

—— = 1 -I- 2ti2 cos {2(j)) (3) 

d(j) 

where (j) the angle between the particle and the 
reaction plane (impact parameter vector), and V 2 
is the elliptic flow. A large V 2 indicates high pres¬ 
sures and early equilibration . Figure ^ shows 
the elliptic flow as a function of centrality, here 
given in terms of Nch/Nmax, where Nch is the 
charged particle multiplicity relative to the max¬ 
imum multiplicity Nmax- The flow is large, and 
is very close to the predictions of hydrodynamic 



Figure 6. Elliptic flow (V 2 ) as a function of pr- 
For pt < 2 GeV/c, hydrodynamic models fit 
the data well. At higher energies, hydrodynamic 
models fail, as expected, and a parton description 
may be more appropriate. The data indicate that 
very high parton (mostly gluon) densities are re¬ 
quired to produce the observed flow. From the 
STAR collaboration 


models that treat the system as a fluid. Flow 
has been studied for identified pions, protons and 
kaons (K^ and ATs), and A. The pr dependence 
of these species flow matches the predictions of 
hydrodynamic models quite well[^. This fluid¬ 
like behavior is another indication of a strongly 
interacting system. 

6. Chemical Fheezeout 

A key observable of chemical freezeout are the 
abundances of various particles. If the system is 
in chemical equilibrium, the abundances should 
scale as exp where m is 

the particle mass, p is the chemical potential of 
the particle (due to it’s baryon and strangeness 
content), k is Boltzmann’s constant, and T is the 
temperature [p8|. Figure ^ compares the ratios of 
a number of different particles, compared with 
the thermal model predictions. The fit finds tem¬ 
perature T = 187 ±8 MeV, baryochemical poten¬ 
tial /ib = 39 ± 4 MeV, strange chemical potential 
Ps = 1.8 ± 1.6 MeV, and strangeness suppression 
factor 7 s = 1.00 ± 0.05. This is hotter than at 
thermal freezeout, which is a later, cooler stage 
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Figure 7. Measured particle multiplicity ratios 
vs. the results of a chemical fit model for the 6% 
most central gold-gold collisions at 130 GeV. The 
model axis errors are the uncertainties in the fit 
output 1^. 


in the evolution of the system. The fj,b is much 
less than the proton mass, showing quantitatively 
that the central region is effectively net baryon 
free. 

The chemical equilibration of hadrons contain¬ 
ing up, down and strange quarks (qg = 1) is 
very different from the situation in e+e“, pp 
and pp collisions. In these elementary collisions, 
7 s « 0.3; strange hadrons are suppressed by a 
factor of 3 below the equilibrium expectations. 

This strangeness equilibration was long-ago 
proposed as a signature of the quark-gluon 
plasma. In a QGP, reactions proceed quickly, and 
equilibrium should be reached rapidly. Reactions 
are much slower in a hadron gas and there are 
many species to produce, so equilibration takes 
much longer. Most (but not all) calculations indi¬ 
cate that the hadron gas equilibration takes much 
longer than the expected system lifetime)^. 

7. Initial States/High Particles 

Study of the evolution of the system before 
chemical freezeout requires a probe particle that 


is created early in the collision. A few probes, 
such as direct photons escape the medium with¬ 
out interacting, and provide information about 
the process that created them. Others, such as 
charmonium and high pT particles, interact with 
the medium and can provide information about 
how it evolves. These probes come from the 
hadronization of high px quarks and gluons (par- 
tons). As of this conference, RHIC has so far only 
presented results on the high px hadrons. 

High px partons are produced very early in the 
collision, at a time r ss li/px- The partons are not 
expected to hadronize until much later, around a 
time t = h/A, where A « 300 MeV is the typical 
QCD scale. Usually, the partons will have exited 
the medium before this point, so that hadroniza¬ 
tion occurs in free space. The medium will inter¬ 
act with the produced parton, not the final state 
hadrons. Since the final state hadron momenta 
depends on the parton momentum, any energy 
loss by the parton in the medium will be reflected 
in the high px hadron spectrum. 

The parton energy loss can be studied by com¬ 
paring hadron momentum spectra from central 
heavy ion collisions with spectra from peripheral 
heavy ion collisions and pp collisions. Published 
results, using the 130 GeV data, have used a pp 
reference spectrum derived from 200 GeV pp col¬ 
lision data from the UAl experiment. 

Fig. HQ, compares charged hadron and 7r° 
production in central gold-gold collisions with a 
normalized pp spectrum. Raa is the cross section 
ratio for gold-gold to pp collisions, divided by the 
number of nucleon-nucleon collisions in the gold. 
In the absence of nuclear effects, Raa = 1- At 
low Px soft (non-perturbative) physics is expected 
to dominate, leading to Raa <?C 1, as observed. 
However, at higher px, where perturbative QCD 
applies well, we expect Raa = 1. In contrast, at 
high Px, in the data, Raa flattens out at about 
0.5 for charged hadrons, and 0.3 for 7r°. 

Nuclear effects, such as shadowing, multiple 
scattering and the net isospin difference can af¬ 
fect Raa- In A A collisions, the incident nucleons 
may undergo soft interactions and acquire some 
initial state px before undergoing a hard inter¬ 
action. Known as the Cronin effect, this initial- 
state Px may increase the final state px and hence 



Au+Au Vs^,|.j= 130 GeV 




- (h"+h)/2 


. ttO 


Pb+Pb(Au) CERN-SPS 

a+a CERN-ISR ./ 



I 


I 

/ . 

-.t I . 


binary scaling 


• * I 




_ 


, 


° 0 2 4 


Pj (GeV/c) 


Figure 8. Comparison of charged hadron and 
7r° pt spectra from gold-gold collisions and 
a pp derived reference. From the PHENIX 
collaboration^ . 


the measured Ra^i- Fig. ^ also shows Raa from 
lead-lead collisions at a center of mass energy of 
17 GeV/nucleon. There, for pT > 2 GeV, Raa 
rises considerably above 1. This is dramatically 
different from RHIC, showing a significant effect 
of the higher energy. In fact, the lower energy 
data shows no energy loss, while the RHIG data 
seems to indicate a very large energy lossQ. 

The systematic errors in normalizing the UAl 
and RHIC data vary with the particle momen¬ 
tum, but are in the 35% range, due to uncertain¬ 
ties in luminosities, cross sections, centrality se¬ 
lection, pseudorapidity distribution, etc. Figure^ 
compares charged hadron spectra from peripheral 
and central gold-gold collisions!^; Ra, the ratio 
of hard particle production in central and periph¬ 
eral AA collisions, per nucleon-nucleon collision, 
is always less than 1. At high pT, Ra ~ 0.3. The 
systematic uncertainties in Ra are about 20%. 

More can be learned about hard interactions 
by considering correlations of high-p^ particles. 
The STAR collaboration has studied the angular 
correlations between particles with pT > 4 GeV/c 
and |r?| < 0.7 (the trigger particle) and a second 
particle withp^ > 2 GeV/c |3^. Figure [T^ shows 



Pt (GeV/c) 


Figure 9. Gomparison of the charged hadron 
Pt spectra from central and peripheral gold- 
gold collisions at 130 GeV. From the STAR 
collaboration!^ . 


the azimuthal correlations, as a function of the 
azimuthal separation A^. There is an enhance¬ 
ment near = 0. This correlation has a similar 
strength and width in pp and AA collisions. 

However, at large separations, A(/) « tt, no cor¬ 
relation is observed in the AA data, while a cor¬ 
relation is seen in the pp data. The back-to-back 
correlations in pp collisions are expected because 
jets are usually produced in back-to-back pairs 
known (although many of the produced jets may 
be outside the experimental acceptance). The 
correlations observed in the pp collisions match 
the theoretical expectations, but the jet pair cor¬ 
relations are absent in the AA data 

The major difference expected between pp and 
A A collisions is anisotropic flow; the solid curve 
in Fig. shows the size of the flow contribution. 
Flow cannot explain the difference between the 
pp and A A curves. 

The suppression of high px particles, presence 
of same-side particle correlations and disappear¬ 
ance of opposite side particle correlations are 
all consistent with a strongly interacting system. 
Only partons produced near the surface of the 
system are able to escape and produce jets. When 
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Figure 10. Correlation function for two high- 
Pt particles. The points are the 5% most cen¬ 
tral (smallest impact parameter) 200 GeV gold- 
gold collisions, while the histogram is 200 GeV pp 
data, both from the STAR detector |Q. 


a parton reaction produces back-to-back partons 
near the system surface, one parton escapes, al¬ 
most unmodified, producing a jet. The other par- 
ton goes in the other direction, into the system, 
where it is absorbed. The large flow at high px 
(Fig. I) supports this picture, showing that even 
energetic particles demonstrate collective effects. 

8. Implications for Air Shower Simulations 

This RHIC data can be used to test air shower 
simulations. As several people at the confer¬ 
ence pointed out, the overall charged multiplici¬ 
ties are considerably lower than most predictions, 
including several popular air shower codes. Mod¬ 
els based on separate hard (described by QCD) 
and soft (phenomenological) models seem to work 
best. 

Gluon saturation models, which can affect the 
depth of maximum shower development, and the 
muon content of showers with energies above 10^^ 
eVQ are disfavored, but not ruled out. In 
these models, the gluon density in heavy ions 
becomes saturated, and low—a: gluons may re¬ 
combine, reducing their abundance. They predict 
that the charged particle multiplicity in AA colli¬ 


sions should be lower than the multiplicity scaled 
from pp collisions; data shows the opposite, with 
the AA multiplicity higher than in simple pp scal¬ 
ing. 

Several other effects are likely to be relevant 
for air showers. Strange particles are copiously 
produced, in chemical equilibrium. This might 
affect the muon content of air showers, compared 
to expectations for pp collisions. 

The presence of non-zero net baryon density 
at mid-rapidity shows that there is a substantial 
amount of baryon stopping, even at very high en¬ 
ergies. In a fixed target frame of reference, these 
baryons carry enormous energy, and so this stop¬ 
ping may have implications for the overall energy 
flow in the collision. 

Several very different analyses show that the 
colliding system interacts very strongly, exhibit¬ 
ing collective behavior that suppresses high px 
particle production. Pure-QGD calculations that 
neglect collective effects may over-predict the 
number of high px particles, and hence the shower 
density far from the core. 

The existing RHIG data is for gold on gold col¬ 
lisions, not the lighter ions and protons found in 
cosmic rays and the atmosphere. Interpolation 
between pp and gold-gold collisions is not easy. 
In the next few years, RHIC will collide lighter 
ions; until then, the various Monte Carlo codes 
can only be tested with light (proton) or heavy 
(gold) systems. 

9. Conclusions 

RHIC is just beginning it’s study of ultra- 
relativistic heavy ion collisions. However, already 
we see a few surprises: large elliptic flow, suppres¬ 
sion of high Px particles, and the complete chem¬ 
ical equilibration of strange particles. Chemi¬ 
cal and thermal equilibrium appear to have been 
reached. This data shows that the system inter¬ 
acts strongly and appears to equilibrate early in 
the collision. 

But, is this the quark-gluon plasma? The evi¬ 
dence suggests that the high densities are strong 
interactions are consistent with a quark-gluon 
plasma. However, a very high density hadron gas 
cannot yet be ruled out. 
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